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Abstract 

The capabilities of one and two-dimensional numerical solar cell modeling programs 
(SCAPlD and SCAP2D) are described. The occasions when a two-dimensional model is 
required are discussed. The application of the models to design, analysis, and prediction 
are presented along with a discussion of problem areas for solar cell modeling. 


Introduction 

Accurate numerical models for single crystal silicon solar ce’ls have proven to be very 
reliable in the simulation of the performance of these cells. These models have proven to be 
extremely useful in: the interpretation of experimental measurements; the identification of 
processes which limit cell performance; the prediction of benefits which will result from 
design and materials changes; the comparison of various cell designs; and the prediction of 
efficiencies which may eventually be obtained in silicon solar cells as various technological 
barriers are overcome. 

The capabilities of a one-dimensional (SCAPlD) and a two-dimensional model 
(SCAP2D) are described and examples of their use for each of the above p rposes are 
given. 1-3 It will be shown that there are circumstances under which cells which appear to 
be one-dimensional require a two-dimensional model to properly simulate their behavior. 

As cells become more efficient the requirements on the accuracy of the physics used in 
the model become more stringent. Effects which are of little significance in poor or 
moderately good cells can take on major significance in high efficiency cells. A number of 
problem areas which are of concern in the modeling of high efficiency cells are discussed. 
These include heavy doping effects, metal-semiconductor boundary conditions, minority car- 
rier mobilities, high injection lifetimes, and carrier-carrier scattering. Each of these may 
have a major impact on the performance of the cell under certain operating conditions. 


The Model 


Physical Effects of Importance 

One of the major advantages of a numerical model is that it affords one the opportun- 
ity to include the very large number of physical effects which may be acting simultaneously 
within a solar cell The complexity of the phenomena and their interactions with each other 
preclude analytic solutions in anything except highly idealized situations, which are not 
indicative of actual cells or operating conditions. An attempt has been made in the formu- 
lation of SCAPlD and SCAP2D to include as many of the physical effects which are known 
to influence cell performance as possible and to do this in a manner which represents our 
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present knowledge of these effects. One of the goa)« in writing these codes was to have 
them be sufficiently accurate in their representations so that they could be used in a predic- 
tive mode. This is possible only if all of the pertinent physical effects are included. 

In those cases where the physics is questionable, we have attempted to include options 
which allow one to choose between various models. For example, in the case of heavy dop- 
ing effects, one is able tc choose between the models of Slotboom, Lanyon-Tuft, and 
Mahan, or to supply a subroutine of one’s own choosing. 

We have attempted to choose materials parameters which in our estimation are the 
most reliable. These materials parameters are used as default values. The user can easily 
change these parameters to values that he views as more reasonable. 

The following physical effects are included in the codes: hole and electron mobilities, 
including ».iieir doping and temperature dependencies; heavy doping effects, using the for- 
mulation of Lundstrom, Schwartz, and Gray; absorption coefficients, including their tem- 
perature dependence; recombination, including Auger, Hall-Shockley-Read, and surface 
recombination. Surface recombination is handled through the specification of the surface 
recombination velocity. In the case of SCAP2D, the effects of surface potentials are also 
included. 

Semiconductor Equations 

The programs perform a full simultaneous numerical solution of the two continuity 
equations and Poisson’s equation subject to the boundary conditions appropriate to one and 
two-dimensional cells. The equations are formulated as shown in equations 1-3. 


V 2y =-i(p- n +N D -N A) , (1) 

V-J p = q(G-R), (2) 

V-J„ = q(R-G). (3) 

The generation term in equations 2 and 3 are given by 

oo 

G(x) = / o (4) 


and the recombination term is given by equations 5. 6 and 7. 
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The hole and electron current densities which appear in equations 2 and 3 are given by 
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J P - -q^pPVVp-kTAipVp 

Jn = -q^ n nVv n + kTAi n Vn 

v p = v-(l-^ (10) 

Aq 

v n=V+7— . (11) 

where v p and v n are the effective potentials defined in equations 10 and 11 and A G and 7 
are parameters which account for variations in the band structure, such as density of states 
ard band gap, and account for Fermi-Dirac statistics. 

No low injection assumptions are made. The equations are solved from contact to con- 
tact with appropriate boundary conditions so that the solutions are valid for all ranges of 
operation and include minority and majority carrier flow. The latter places s°me restric- 
tions on the CPU word size required for solution. 

These codes have been extensively tested for accuracy by comparing the results of 
their predictions with experimental results obtained on very carefully and extensively 
characterized cells for a wide range of cell designs and operating conditions. The agreement 
has been such that a high degree of confidence has been developed in results 'mputed 
using these codes. 


Code Description 

Figure 1 is a block diagram of the structure of SCAPlD and SCAP2D. The operator 
must supply information about the materials parameters, a description of the device to be 
analyzed, the type of analysis which he w : shes to perform, and the spectrum, if appropriate. 
He also can, if he wishes, control some of the details of the numerical solution; the amount 
of information supplied while the program is converging to an answer and how the output 
information will be stored or displayed. 

The results of the computation are presented in printed summary form and the 
detailed results of the calculation are stored on magnetic tape. A separate plotting routine 
is used to access the information on tape and to display the appropriate parameters. The 
plotting capability is one of the most valuable features of the code, in that it allows one to 
effectively have a microscopic view of most of the parameters of interest in the interior of 
the cell under operating conditions. We will show some of the available graphical output as 
we discuss the capabilities of the code. Table I shows the input control offered to the 
operator. In every case default parameters are specified if the operator chooses not to sup- 
ply a parameter. 

Table II contains a listing of plots which are available through the plott ng program. 
In this case the operator specifies the type of plot which is required and the region of the 
cell for which he desires that plot. Most of the figures which follow were obtained directly 
from this plotting routine. 

In addition to the reliability of the output, the utility of codes of this type will depend 
on their ease of use and efficiency of computation. For example, in a design mode, it is 
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Figure 1 Block Diagram of the Structure of SCAPlD and SCAP2D 


Table I - Input Parameters 


Device Description 
Doping Profiles 
Step junction 
Erfc (N s ,Xj) 
Experimental Profile 
SUPREM n 
Dimensions 

Materials Parameters 
Lifetime (r and energy) 
Surface Recombination 
Auger 

Bandgap narrowing 
Slotboom 
Mahan 
Lanyon-Tuft 
User supplied 


Spectral Choices 
AM 0 
AM 1.0 

AM 1.5 direct & global 
Monochromatic 
Uniform generation 
User supplied 

Optical Filters & Reflection 
Filter (Ge, Si, Si0 2 , GuAs) 
Back surface reflector 

Types of Analysis 
Dark I-V 
Illuminated I-V 
Solar Cell 
Spectral Response 
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Table II - Plotting Options 


Carrier concentration 

Hole and electron current -nsities 

Change in potential (from qu ibrium) 

Doping density 

Energy band diagram 

Electric field 

Hole an - electron quasi-electric fields 
Effective fields (electric plus quasi-electric) 
for holes and electrons 
Optical generation 


Hole current density and components 

Electron current density and components 

Mobility 

Lifetime 

Ratio of n lf /n io 

Potential 

Recombination rate 

C harge density 

Excess carrier concentration 


advantageous to be able <o make multiple runs in a reasonable length of time and at rea- 
sonable cost. Whiie S CAP ID can be run effe lively on nearly any mainframe computer (a 
typical run on a CDC 6600 require* 100-300 CPI seconds). SCAP2D requires a very fast 
machine with a large amount of actual or virtual memory. On a Cyber 205, 300 CPU 
seconds are required for a typical run. 


Situations Requiring Two-Dimensional Analysis 

In man\ situations a one-dimensional simulation is quite adequate and there is no need 
to use the more complex and exnensive two-dimensional simulation. On the other hand, 
there are a number of situations which only a iwo-dimensional simulation will suffice. 

Some of the situations which require two-dimensional analysis are quite obviour, while 
others appear a be one-dimensional in nature, hut. in fact, require a two-dimensional solu- 
tion for proper description of the cell performance. Most of the cell structures which have 
been proposed as high efficiency silicon cells fall into the obviously two-dimensional analysis 
category. Among these structures are the Interdigitated Back Contact cell, the Vertical 
Multi-Junction cell, the Etched Multiple Vertical Junction cell, the Polka Dot cell, and the 
Grating cell. .As an example of the use of SCAP2D in the analysis of these two-dimensional 
cells, we show figures 2 through 4 for an IBC cell. In Figure 2 we show the total short cir- 
cuit current flow under one sun-conditions. In Figures 3 and •! we show the majority and 
minority carrier flows for this same cell operating under the same condition. 

Less obvious applications of the two-dimensional code are shown in Figures 5 through 
6 in which a conventional solar cell has been analyzed. In Figure 5 we show the potential 
distribution along the emitter from a poi:.! half way between the grid lines up to the grid 
lines under open circuit conditions. This figure illustrates that there is a lateral voltage 
drop along the emitter, even under open circuit conditions, as a result of the current which 
is injected in the vicinity of the grid line. Figure 6 shows the circulating currents which 
exists in the vicinity of the grid ii *\ 
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Figure 4 Minority Carrier Current Flow for the Cell of Figure 2 
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Figure 5 Potential Distribution in the Emitter of a Conventional Solar Cell Operating at 
•100 suns. 
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Figure 6 Circulating Current in the Vicinity of a Grid Line for a Conventional Solar Cell 


An even less obvious problem with on e-dime ns ion -al simulation occurs when one tries to 
properly model the front surface of a conventional cell. This surface is composed of a 
SiOj-Si interface and a metal-semiconductor contact. In a one-dimeDsionaJ simulation, one 
is forced to aggregate the two effects with some equivalent front surface recombination velo- 
city, Sp Figure 7 illustrates the difficulty with this approach. Under short-circuit condi- 
tions the proper value of S F is equal to the surface recombination velocity of the 
interface. However, near op*>n circuit cond'tions., the proper v of Sp may be 3 to 4 ord- 
ers of magnitude larger. This is a result of the fact that the ..! semiconductor contact 

may be a very effective recombination site for minority carriers. . is particularly important 
as the operating voltage of the cell increases. For proper operation of a one- dimensional 
code, the front surface recombination velocity ahou'd be a function of operating condition. 
The two-dimensional code does not have this problem, since the surface recombination velo- 
city at the SiO^-Si interface and the metal semiconductor interface are specified separately, 
and the recomuination along the entire surface is properly accounted for under all operating 
conditions. 

At high operating conditions, such as are found in concentrator solar cells, even the 
conventional ceil behaves in a two-dimensional fashion and must be modeled using the 
two-dimensional code. Minority carrier current flow for a conventional cell operating at 800 
suns is shown in Figure 8. If this cell is modeled using the one-dimensional code under 
these operating conditions, serious errors are encountered in the computation of the fill fac- 
tor which can not be compensated for by including an external series resistance in the 
model, as the effect is nonlinear. 
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Figure 7 Effective Surface Recombination Velocity as a Function of Operating Voltage and 
the SiOj-Si Surface Recombination Velocity 




Figured Minority Carrier Current Flow for^a Conventional Solar Cell Operatir at 800 
Suns (V = .600 volts J =21.6 ainp/cm 2 ) 
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Modes of Utilisation 


As we mentioned previously, a carefully prepared solar cell model is useful in a number 
of modes. In this section we will discuss the use of SCAPlD and SCAP2D as a design tool, 
a sensitivity analysis tool, an aid in the analysis of experimental data, an aid in the provi- 
sion of insight into the operation of the cell, and, finally, as a predictive tool for the com- 
parison of proposes cell designs and as a means of projecting performance as various tech- 
nological barriers are removed. For the sake of continuity, we have chosen to use the San- 
dia nigh concentration cell operating at 1 sun as a base line design. This is a cell which has 
exhibited 18?c conversion efficiency ?t one sun, and 20^ conversion efficiency in the 50-100 
sun range for an AM 1.0 spectrum. 


Design 


As a simple example we show, in Table III, the effects of variations in the base doping 
about the present design doping of 2.29 x 10 16 , on the performance of this cell. We see that 
the present base doping is nearly optimum for the design parameters used in the other 
parts of the cell. 


Table m 

Solar Cell Performance Dependence on Base Doping 
AM 1.0 (one sun) 


Base Doping 
cm - * 

V« 

volts 

I sc 

ma/cm 2 

F.F. 

Efficiency 

% 

5 x I0 15 

.634 

35.1 

.828 

18.35 

1 x 10 16 

.640 

34.8 

.833 

18.46 

2.29 x 10 1 * 

.649 

34.4 

.836 

IS.55 

1 x »0 17 

.656 

33.3 

.838 

18.21 

5 x 10 17 

.650 

30.2 

.836 

1637 


Sensitivity Analysis 


By utilizing a computer code such as SUPREM to simulate fabrication conditions one 
can model the sensitivity of device performance to fabrication parameters. Here, as an 
extreme case, we examine the effects of changes in the emitter doping profile on cell perfor- 
mance. The Sandia cell was simulated using the two emitter profiles shown in Figure 9. In 
Table IV, a comparisoa of these simulations is shown. Note that the erfc emitter profile 
simulation predicts a higher V oc . This is due to the lower net recombination in the emitter 
as compared to the SUPREM II emitter profile simulation, as shown in Figure 10. Recom- 
bination is higher in the SUPREM II emitter because the doping is higher over most of the 
emitter volume, and therefore Auger recombination is correspondingly higher also. 

If the results of a process simulation program such as SUPREM are coupled with 
SCAPlD or SCAP2D as shown above the sensitivity of the cell to process variations can be 
readily established. 
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Figure 9 Emitter Doping Profiles as Determined by SUPREM D and Complimentary Error 
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Figure 10 Emitter Recombination for the Two Doping Profiles Shown in Figure 9 
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Table IV 


Dependence of cell performance on emitter doping profile 
AM 1.0 (one sun) 

Type of Profile 

volts 

V * 

rna/cm 

FF 

Efficiency 

% 

Erfc 

SUPREM n 

.648 

.632 

34.3 

33.0 

.836 

.833 

18.6 

17.75 


Analysis 

It is possible, by adjusting the parameters entered into the code, to obtain a fit 
between the model generated results and experimental dark I-V, solar cell, and spectral 
response curves. If this fit can be obtained for a single set of parameters, then one has a 
reasonable expectation that these are the correct parameters describing this device. 

Insight 

With the ability to observe most of the parameters of interest as a function c* position 
and operating » ditions anywhere within the cell, it is possible to achieve a gru:, deal of 
insight into the uniting factors on any cell design. Examination of the model for. the 20 % 
Sandia cell very quickly establishes that the cell appears to be emitter li.nited, and, in fa. . , 
that further efforts in improving the performance of the cell should be devoted to reduction 
of the metal-semi condur tor contact recombination and in reducing the volume of the 
heavily doped emitter. 

Prediction 

Potentially one of the most valuable, and also one of the most risky, uses of the 
numerical models is as a predictive tool. The models have already been shown to be quite 
reliable in comparing the relative merits of different cell designs. One particularly attrac- 
tive way to utilize the code is to use it to identify limiting phenomenon in a particular cell 
design and then to remove that limitation and observe the effect on cell performance. In 
this fashion, one can predict benefits which will accrue through various advances in technol- 
ogy, and, in fact, can make some reasonable estimates of the ultimate performance of sili- 
con single crystal solar cells. This latter use of the code is particularly risky since as the 
performance of the cell improves, physical effects which may have been insignificant in their 
effect on cell performance before, may suddenly become the dominant limitation. 


Problem Areas 

There are a number of areas in which there is concern about existing solar cell models 
either because the physics is not well understood, available data is thought to be unreliable, 
or because the effect has not been include in the model. These areas of concern are dis- 
cussed below. 



Heavy Doping Effects 


There is a controversy over the origins and magnitude of heavy doping effects. There 
is a great deal of scatter in the measured effective band gap narrowing, particularly in the 
very heavily doped samples where we have our major concern. In order to alleviate this 
situation somewhat, we have provided the operator with the option to choose between most 
of the popular band gap narrowing models. This remains an area of major concern and is 
probably the least reliable area in the modeling of silicon solar cells. 

Auger Recombination 

Some uncertainty exists about the reliability of published Auger coefficients. At least 
two groups (Sandia and General Electric) have indicated that published Auger coefficients 
may be too large. 

Minority Carrier Mobility 

Reliable measurements of minority carrier mobility do not exist. Various authors have 
proposed that the minority carrier mobility is larger, smaller, and the same as the majority 
carrier mobilities of the same type carrier. As a consequence, SCAPlD and SCAP2D 
assume that the minority carrier mobilities for electrons are the same as they would be if 
electrons were majority carriers. A similar assumption is made for holes. 

Metal- Semiconductor Contacts 

In well designed high efficiency solar cells, the metal semiconductor contact limits the 
open-circuit voltage. The lemoval of this high dark current source, through the use of tun- 
neling contacts or through the reduction of the metal-semiconductor contact area, has 
already demonstrated a significant improvement in open-circuit voltage. F urther advances 
in this area may well employ heterojunction structures in addition to the present tunneling 
structures. SCAPlD and SCAP2D allow for specification of a finite minority carrier surface 
recombination velocity to model this effect. 

Doping Profiles 

We have already seen that device pf rfovmance can be a strong function of the shape of 
the emitter doping profile. SCAPlD and SCAP2D allow for the use of a complimentary 
error function, a computed profile based on the Fair diffusion model for phosphorus, doping 
profiles obtained from a process simulation program such as SUPREM, or experimental 
data. The use of data from SIMS measurements has the problem that it includes the total 
impurity concentration not just the electrically active dopants. If any precipitation is 
present in the highly doped region, SIMS will overestimate the amount of active dopant. 
Spreading resistance measurements are a measure of the free carrier concentration. Near 
the depletion region this can lead to significant errors in the doping profile if the spreading 
resistance profile is interpretated as being the same as the doping profile. 
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Spectral Response 

Spectral response measurements are particularly useful for obtaining information about 
the base lifetime and the surface recombination velocity. However, some difficulty is 
experienced in matching long wavelength response with computed response. This difficulty 
has been traced to the fact that small changes in device temperature can lead to large 
changes in the long wavelength response as a result of changes in the absorption coefficients 
due to a shift in the band edge. 

In order the determine the surface recombination velocity cf the Si0 2 ~Si interface, it is 
also highly desirable to have spectral response measurements in the very high absorption 
regime of .35 - .4 jim. Accurate measurements of the internal quantum efficiencies are 
difficult to obtain at these wavelengths. 

Effects of Band Gap Narrowing on Long Wave Length Absorption Coefficients 

At the present time no corrections for the effect of band gap narrowing are made to 
the absorption coefficients. 

Carrier-Carrier Scattering 

Carrier-carrier scattering can be a significant effect in high concentration solar cells, 
and will become a significant effect in one sun solar cells as the efficiency is increased. 


High Injection Lifetime 

At the present time very little data is available on majority carrier lifetime. A typical 
modeling approach is to assume that the majority carrier lifetime is the same as minority 
carrier lifetime. This seems to give reasonably good agreement with cell performance under 
high injection conditions, but direct measurement of the high injection lifetime would be 
highly desirable. 


Conclusions 

One and two-dimensional device models have been quite successfully employed as an 
aid to design, interpretation, sensitivity analysis, and prediction. However, the predictive 
capability of any device code is only as good as the physics which is modeled and the data 
which is supplied. If further improvements are to be made in the performance of single cry- 
stal silicon solar cells, careful attention will have to be paid to both of these areas and & 
great deal of effort will have to be devoted to measurement techniques which will allow the 
independent determination of the parameters which must be supplied to the device code. 
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DISCUSSION 


SAH: I would like to ask you about that particular example you gave, the com- 

plementary error function, also the SUP2EM. There is quite a difference 
in the results. What is the basic mechanism? What are the recombination 
mechanisms t’jat give you the difference? 

SCHWARTZ: What's included in the code is Hall-Shockley-Kead and Auger using 

Schmidt Auger coefficients, and for that particular run — thank 70U for 
asking — 1 should have mentioned at the time: for that particular run we 
used the Slot-Boom band-gap narrowing model, so there is a significant 
amount of band-gap narrowing occurring, and recombination echanisms. For 
that particular run. Auger and Hall-Shockley-Read, and I didn't bring the 
plot along. The plot shows the split between them. I don't recall what 
the split was. 

SAH: For the particular profile, was it the Auger that causes the one to be 

better than the other one? 

SCHWARTZ: Yes, 1 believe that it was Auger, but I don't have the plot with 

me. 


NEUGROSCHBL: You said that the published Auger coefficients don’t agree with 

the experiment. In order to get agreement, do you need larger lifetimes 
or shorter lifetimes? 

SCHWARTZ: The recombination wants to be reduced. I should qualify that a 

little bit. I was repeating what is said in a couple of publications by 
Posene at GE and Weaver at Sandia. It is possible that the problem lies 
in the band-gap narrowing model and not in Auger, so one wants to be a 
little careful. There is a problem in the emitter, and that's clear, and 
people have tended to blame Auger. I guess I'm not completely convinced 
that it is Auger — it may <*ell be related more to the band-gap narrowing. 

DAUC: I would like to follow up on Sah's question. You have the same carrier 

density at the surface for the SUPSEM and for error function. So, norm- 
ally you have much larger field right ac the surface in case of comple- 
mentary function, and I would expect less recombination there. Would you 
give some reason why? 

SCHWARTZ: Yes. There are a number of reasons. One is that the recombination 

v only depends on the Auger coefficients and lifetime, but also on the 
excess minority carrier concentration, and if I had shown the plot, what 
you would have seen is that aany of the carriers are recombining in the 
case of SUPREM as they moved, and the axis cure was lower at the surface. 
The other difference is in the way that the band-gap narrowing effective 
field is distributed. The minority carriers in the emitter don't see 
just an electrostatic driving force due to the gradient. There is another 
component, which is associated with band-gap narrowing itself, and it 
tends to reduce the effect of pulling minority carriers away from that 
surface or keeping them out of the emitter. Both are operating and they 
are distributed differently. 
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DAUD: The second question has to do with actually running this program. We 

find that one of the items that one has to put in is the tpQ and 
T n o* Normally when we measure, we either measure the lifetime or 
diffusion length where the doping is already there. How does one 
reconcile with this? What kind of numbers one should put for t_q 

SCHWARTZ: I'm sorry. You say you noraally measure 'That? 

DAUD: Say we hav~- a cell where we measure the diffusion length or the life- 

time. We cannot directly put that in your model because y^ur model 
corrects it for the doping. 

SCHWARTZ: Yes. There were a lot of features that were not talked about here, 
and since we have sent a copy of the code out to JPL, he is asking. You 
have the option of turning on what amounts to a Kandel fit to doping. If 
in fact you have that turned on and if the base doping is above the 
transition doping fot the the Kandel fit, then you have to correct that. 

If you have measured the lifetime at that particular doping level, you 
either have to correct it or just turn Kandel fit off and enter the value 
you measure. It's your choice. It is under your control on the input 
deck. 

LIMDKOLM: I have three questions. They are all, I think, fairly quick. Just 

to remind the audience: you recast some of the -- what you might say, 
equations which were truer to the physics — into a form that one is more 
used to seeing in a conventional treatarent of semiconductor device 
physics. In so doing, you introduce the parameter that you called small 
gasma, lower case gamma, and I think that that parameter was supposed to 
have taken care of various things that were being violated by the density 
of states. You made a big point, which is an extremely valid point — 
that the parameters that go into the model have to be measurable. So 
with that preparation for the audience, I could have just asked you what 
success havw you had in measuring gamma, and how do you do it? How don't 
you do it, if you can't do it? 

SCHWARTZ: I tried to stay away from the equations, Fred, and I apologize for 

the poor quality here. The gamma that Fred is talking about entered into 
the effective potential that we showed. Here is the electrostatic poten- 
tial, and here is the term iiich I said was an effective ayrraetry factor, 
and that term is not normally known. As one gets measurement of the 
effective band-gap narrowing, which in fact looks like this — the band 
gap plus all the degeneracy and band structure effects. Gamma, in fact, 
has electron affinity divided by the delta G minus terms for degeneracy. 
The answer to your question is, you don't know. But before I let that 
go, it turns out that for solar cells you seldom care, and the reason for 
that is the following: GB observed this first, the range on that nimbers 
from 0 to 1 and one can run the full range and see almost no deteu jble 
change in the device characteristics. On the other hand, if you look 
inside the device, there are radical differences in the electric field 
distribution ?n the emitter, in that region. There are huge differences. 
But it turns out, and you can do this in closed form, that if you are 
dealing with a region which is quasi- neutral and low injection, as the 
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emitter is under anything except extreme conditions, you can show in 
analytic form that, in fact, terminal characteristic is independent of 
tbit parameter. There are places where it makes a difference. It does 
make a difference in junction capacitance: you can show that fairly 
readily. It makes a difference in the electric field. And a dramatic 
one if, for instance, you ar' concerned with avalanche, not in a solar 
cell, but calculations of avalanche, and you are in a heavily doped 
condition. Then possibly one should be worried about present-day 
avalanche coefficients if they are based on classically calculated 
electric fields. But the answer is that it doesn't affect the terminal 
characteristics. 

LINDHOLM: The reason I asked that is twofold. First of all it would seem to 

me that for diagnostic purposes it would be desirable to k^ow what you 
call gamma, or what I would rather just call the electro. finity, and I 
know about the comment you made about tie quasi-neutral ’. p ion. But as 
you start entering a little bit into the junction transition region, then 
I think it becomes more important. And the reason 1 asked that is, if I 
correctly read your earlier paper with Mark Lundstroji, that you indicated 
a method for measuring gamma. And so I come back to mv original question: 
what degree of success have you had in measuring gam^a? 

SCHWARTZ: Very little. We're still working on it. 

LINDHOLM: T think it is a good thing to work on, actually. The reason is that 

it's very easj, even though you did not intend to do this — in fact, your 
wording was very careful — but people will take sort of '.ial cases and 
say it doesn't matter, but it does matter from a diagnostic point of view 
in finding out what's going on with the profiles ,how you can improve the 
device, that kind of thing. I'm very glad you used such careful wording, 
so congratulations on that. 

SCHWARTZ : I'm glad you read the paper, v „ry few people have. 

LINDHOLM: Its a very interesting and very good paper. I was extremely inter- 

ested in your measurement of gamma. 

SCHWARTZ: I have a Ph.D. student who is extremely inter sted too. 

LINDHOLM: The other thing that I noted that you said, an:* since you word 

things so carefully, I was noticing that you said that most of the people 
who made electrical measurements in effect were measuring the p-n prod- 
uct. Now, I think that that's probably true of the Slot-Boom graph you 
use, the transistor structure; I think it's not true of ERIC people 
Would you agree with that? A good portion of the data now c^.t’.ig out is 
EBIC data and the guys from GE who aren't here **re talking a y p-n 
product and I think that they can't do that. The fundamental reason is 
that Fermi levels have to remain sensibly, spatially 4 ovarian.: ?’ jr a 
significant region of the device in order for that measurement tc yield 
the p-n product. 

SCHWARTZ: That's absolutely right, and from a physics point of view is ver: 

pleasing. And from a modeling point of view it's difficult, because now 
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in the intermediate region, where one worries about what happens to the 
band shape, yon * ve got a problem in using it. 

LTBDHOLM: One last one. You keep mentioning two-dimensional programs, and I 

would like to know about three-dimensional programs, and why you can get 
away with two-dimensiona 1 . 

SCHWARTZ: It’s not so much getting away with two-dimensional, it's really a 

level of difficulty. In our r»se -- one thing £ didn't mention, and I 
neglected it inadvertently -- i. told you that it cost about $5.00 and 
takes 100 to 300 seconds on a 6500 CDC if you run the 1-D code. It takes 
about 300 CPU/second on our CYBBR 205 Supercomputer to run the 2-D code. 
Mow, what I mean by that is to run a full analysis, something like 10 
points, and do all the associated stuff. That uses about 2,000 mesh 
points. So right now it's about what we are capable <_f tackling. If 
Purdue will put another two million words into the main nemory of that 
computer than we will look at the three dimensional. 

LIKDHOLM: If the next speaker is successful in cutting the computer time 

. iguif icantly by his technique, as he suggests, would you then strongly 
advocate three-dimensional steady-state modeling as a highly useful, 
moderately useful rarely useful vehicle for solar cell design optimi- 
zation and for mot : -oring manufacturing processes? 

SCHWARTZ: I think that. * very good question, whether it's asked about one, 

two or three-dimensional, ana the answ"* 1 Les in how easy it is to use, 
whether it's fast, are the tui*-- arounds quick, *nd what is reasonably 
''heap. Because you do have to make a l^t of runs, and if it's very 
expensive or very time- consuming, the utility becomes a lot less. 

LIMDttOLH: Suppose it doesn't cost anything? 

SCh .ARTZ: Then its very useful to do even for one sun. Is that what you 

think? 

LIMDttOLH: I haven't thought enough (.bout it. 

SCHWARTZ. If you want to do a cell like Dick Swansou e. 

LIND HOLM: A more conventional cell is very useful there? Or is it moderately 

uceful' 

SCHWARTZ: I doubt it, I can't see the benefit to a conventional cell with 

thvee dimensions. 

QUESTION: Dick, a quick question. Did you decide that radiative recombina- 

tion and trap Auger effect could be neglected safely? 

SCHWARTZ: bo, we didn't. I told you the status of the code as it was. If is 

r fairly straightforward matter to add those components to it, we just 
haven't driven things, we haven't mad'' any _uns where the other lifetime 
limiting met 3*s were low enough to dc that. But clear.j , if you do, 
that s a limi.iug mechanism that is not present and needs to be added. 
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QUESTION: And op the surface recombination you might a 1 so have trap Auger 

effects, or do you wrap it all up in a surface recombination velocity? 

SCHWARTZ: Well in the one-dimensiouel code it's wrapped up in the surface re- 

combination velocity. I didn't talk about how it's handled in two 
dimensions. In two dimensions you eould either do it by integrating 
through the trap states and capture cross-sections, or you could do it by 
a lumped parameter, which is probably not as good -- and you do have to 
control surface potential, which we do by setting the charge in the oxide. 

QUESTION: Now, a last question: capacitance. Do you work it out or do you do 

current- volt age, capacitance voltage? 

SCHWARTZ: Mo, JPL doesn't know if their version of the code does have capac- 

itancv in it, we just didn't tell them. It is the equivalent of very 
low-frequency capacitance I'll tell you what it is and you can name 
it. It is the integral of either the electrons or holes with the appro- 
priate voltage term put in, and you are quite right. 

QUESTION: It can miss by a factor of three or four? 

SCHWARTZ: Yes. We don't use it that way. 
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